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Abstract 

We analyze the /3-functions of Yukawa and electromagnetic theories with Lorentz vio¬ 
lation (LV) and propose an alternative method to find the scale dependence of the different 
fields that parameterize such violations. The method of solution consists of decomposing 
a family of parameters into their irreducible representations and thus generating a group 
of subfamilies that obey the same symmetries and transformation rules. The convenience 
of this method relies on the fact that the degree of complexity of the problem is reduced 
in the decoupling process. For a set of parameters describing a Lorentz-violating the¬ 
ory, we expect their associated /3-functions to be nonnegative or, otherwise, their scale 
dependence weak enough. These conditions are necessary because asymptotically-free 
parameters could leave high imprints of LV at low energies, which are ruled out by ob¬ 
servations. Besides imposing some constrains on the coefficients that describe LV, this 
method can be used to fold the fixed points of the theory in each particular irreducible 
subspace. 


: a.ferreroOuniandes.edu.co 



1 Introduction 


Although LV has not been observed in our relatively low-energy experiments, there is no 
reason to assume that small violations of Lorentz and CPT symmetries are incompatible 
with quantum field theories met Nonetheless, any model that includes Lorentz or CPT 
violations, for instance, the Standard Model Extension (SME) [HE], [3], must be compatible 
with the phenomena observed in our universe. 

The conservation of Lorentz and CPT symmetries should be related. In 2002, Green¬ 
berg proposed that CPT violations imply the violation of Lorentz invariance [I]. Never¬ 
theless, Deutsch and Gracia-Bondia recently suggested that such claim is still on some¬ 
what shaky ground |5j. While the relation between both symmetries might still be 
unclear, there are some fundamental causes that could induce their respective viola¬ 
tions. Some possible scenarios include models with time-variating coupling constants 
0 0 El El EDI El] and models with noncommutative geometry, spacetime discretization, 
among others pH E31 EH E5] • 

The values of some of the parameters described by the SME can be constrained by 
means of experiments. Most of the analysis in the SME has been performed in the 
electromagnetic sector [T6| |T?, 18} EH, 20j. Additionally, some work has been performed 
in the gravitational [21, 22] and neutrino regimes [23] • A possible violation of Lorentz 
symmetry in the Yukawa sector has not been strongly emphasized in spite of its great 
importance [24] ]25], specially with the great evidence supporting the existence of the Higgs 
particle coming from the LHC experiments [27]. The existence of the Higgs particle, which 
generates the masses of the particles in the standard model, opens the window to have a 
better understanding of Yukawa interactions. 

Besides using observations to test the validity of models that include Lorentz and 
CPT violations, the renormalization group and the scale dependence of the parameters 
that describe such violations are important to formulate a consistent theory. As it is well 
known, the scale-dependence of a physical quantity in a quantum field theory is described 
by its /3-function. While a positive value of such function implies that the observed value 
of its associated physical quantity increases with the energy scale and a negative value 
describes the opposite behavior, the zeroes of the /3-function provide information of the 
fixed points, in which the physical quantities have no scale dependence. 

The electron charge in QED is a very famous case in which a positive /3-function is 
obtained. In this case, therefore, we expect its observed value to decrease in the low energy 
regime. Similarly, it is natural to expect the same behavior for a family of parameters 
describing LV in the SME. If some of their associated /3-functions turn out to be negative, 
we should expect to observe (unless their scale dependence is weak enough) a relatively 
large imprint of LV at low energies, which, at the best of our knowledge, has not taken 
place. Therefore, we can impose constraints on some of the coefficients contained in 
Lorentz-violating theories by studying their scale dependence; the validity of such models 
can be tested with this analysis as well. 
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Indeed, if the coefficients describing LV turn out to have positive /3-functions, we 
could say that Lorentz and CPT symmetries are low-energy emergent symmetries present 
in our universe |28, [293- While (if the described scenario is correct) the universe and 
its spacetime structure might be described by very different physics and other symmetry 
principles at high enough energies, there could exist an energy scale from which such 
effects are practically washed out. This can help to explain why Lorentz symmetry seems 
to be exact. 

The purpose of this paper is to provide a method to easily solve the /3-functions in 
some theories with LV and impose constrains on such models by demanding the con¬ 
dition that their /3-functions must be nonnegative—this condition can be used to test 
the consistency of such models. Particularly, we will apply this method to Yukawa [30] 
and electromagnetic theories [31, [32] at one-loop order. The /3-functions for a Yukawa 
theory with Lorentz and CPT violations were already found in [30]; however, they were 
not solved. On the other hand, the one-loop order renormalization of Lorentz-violating 
electrodynamics in flat spacetime was already studied in [3T]; nonetheless, irreducible 
representations were not used to decouple the sets of differential equations, allowing us to 
evaluate the convenience of our method. In addition, by separating the Lorentz-violating 
operators into different symmetries, the fixed points of the theories can be studied in 
each particular subspace. It could also be interesting to apply this method to electroweak 
interactions [22], non-Abelian gauge theories [32, 23]. or other scalar-fermion interactions 
[36] with LV; however, this will not be done here. 

This article is organized as follows: in section [2] we briefly describe the Lagrange density 
for a Yukawa theory with Lorentz and CPT violations as well as the different coefficients 
that describe such theory. Section [3] describes the general method of solution used to 
solve for the running couplings. Particularly, such method is applied to Yukawa theories 
in section [2] the application of such method in a known model such as electromagnetism is 
studied in section[5]for comparison reasons. Finally, a brief summary and some conclusions 
are presented in section |6] 

2 Yukawa interactions with LV 

As stated in [30], a theory of Yukawa interactions with Lorentz and CPT violations can 
be described by the Lagrange density 

C Y = \d tl <i>d tJ <i> + - i/i 2 </> 2 - ^A</> 4 + — M) i/i — (1) 
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— 
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O 
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O 

— 

— 

+ 

+ 

— 

— 

+ 

Hij) Lij , djOi doj 

— 

+ 

— 

— 

+ 

— 

+ 

4 

bo, Jo, 9m, Qijk 

+ 

— 

+ 

— 

+ 

— 

— 

bj, Jj, 9i0k, 9ij0 

+ 

+ 

— 

+ 

— 

— 
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Table 1: Discrete symmetry properties that describe the mixing of operators. 


where irrelevant total-derivative terms have been ignored and 

r = 757 m+^+ *75r + , (2a) 

M = m + i^ b m' + 757 ^ + , (2b) 

G = g + * 75 // + + 2 ^ U(T ^ u ' (2c) 

When quantum corrections are introduced, the previous operators mix according to their 
symmetry properties. This fact allows us to introduce a set of groups whose elements 
obey the same symmetries, which are described in table [0 The set of counter-terms and 
the renormalization of the theory is studied in detail in [ 30] . 


3 Renormalization 

In our previous paper [30], we found expressions for the /3-functions associated with all the 
parameters that describe LV in a Yukawa theory. Nonetheless, solving for the running 
couplings was not performed due to the complexity of the /3-functions. Such difficulty 
arises when different parameters mix due to quantum corrections and so generating sets 
of coupled differential equations. 

In our analysis, each group of parameters will in turn be decomposed into subgroups 
using the symmetric group; the action if this group is the permutation of their indices. For 
example, the coefficient dJ w can be decomposed into three irreducible representations: a 
trace-like part written as r/ /w d“ Q , an antisymmetric part described by d'^ = 7j(d*“' — d Ufl ), 
and a traceless symmetric part given by d (Y = |(<3 /iI/ + d ul1 ) — \rj^d a a . 
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It is well know that, in this case, the trace-like representation has scalar (rank-zero 
tensor) properties; the antisymmetric and symmetric parts have vector (rank-one tensor) 
and rank-two tensor properties respectively. This situation also arises, for instance, when 
we add the angular momentum of two spin-one particles. 

If the mixing of the different coefficients in a Lorentz-violating theory did not exist, 
we would obtain /3-functions of the form 

(Ar)' 11 '"'*" = f x {m , rri, g, g')x^ ^ n , (3) 

where f x {m,m !, g, g') is a function that depends on the masses and the couplings of the 
Lorentz invariant theory (notice that the sign of (/3 x ) IJ,1 "' IJ,n depends on the sign of f x ). 
We use the same notation as in [30) [37] by defining the dimensionless quantity p = p/M, 
with M a particular mass scale, to represent the scale associated with the renormalization 
group. The variable or x^ 1 '" ,, ' r ' (p) will indicate the value of such parameter at the 

scale p and its value at M (p = 1). For parameters satisfying Eq. ([3]). its scale 

dependence is given by [30] 



Nonetheless, when we attempt to obtain the scale dependence of the Lorentz-violating 
coefficients belonging to the same group, we will obtain coupled systems of differential 
equations of the form 


N 




m 




k =1 


(5) 


where M mk are matrix elements that depend on some functions ( f x )mk and N is the 
number of elements belonging to the same group or subgroup. Each subgroup is in fact 
described by means of a decomposition into irreducible representations, where all elements 
belong to the same irreducible class. The symmetric group S n will be used to accomplish 
this goal. 

If we introduce the diagonal matrix given by M — AM^ A~ l we can decouple 
Eq. (J5J) and obtain solutions of the form 


(d)fli ■■■[In -(d)pi--fin q 

x k ~ x k J kk i 


where 


S kk = exp 



Um,rri,g,g')(p') 

p' 


( 6 ) 

(7) 


are the solutions to the /3-functions in the diagonal basis. Here are the eigenvalues 
of the matrix M; the initial conditions in the diagonal basis are related to those in the 
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original basis by = A 1 x /il ”' Mn . Therefore, the solutions, in the original basis, 

are given by 

x w-^(p) = ASA- 1 ^ 1 '"^ . (8) 


Notice that the set of eigenvalues determine the sign of the /5-functions in the diagonal 
basis. However, when we have both positive and negative eigenvalues, the coefficients in 
the original basis have contributions that grow and decrease with the energy scale. This 
still predicts a large imprint of LV at low-energy scales unless some very special conditions 
are satisfied or the scale dependence is weak enough. In the presence of fixed points we 
have null eigenvalues, which indicates that the /5-function of at least one parameter can 
be written as the linear combination of the other ones. 

It should be emphasized that a decomposition into irreducible representations is not 
the only method that can be used to decouple the sets of differential equations. However, 
this method greatly reduces the complexity of the problem. For example, let us suppose 
that we have the set of /5-functions 

(Pur = \uirTU a a + \v^ v V\ + u u U^ + u l2 U + v n V^ + v 12 V^ , (9a) 

(/5 vY u = \u 2 rTU a a + \v 2 rTV a a + u 21 U^ + u 22 U ^ + v 21 V^ + v 22 V^ , (9b) 

where iq, Uij, ly, and are functions of p. If we permute the indices {y, u} we obtain 
the linear system /3 X = Mf, where x 1 = (U^, U uti , V J ‘ V , V u, \ rj lxiy U a a , rj ,w V a 0 ) and 


Uu 

U\2 

vu 

V\2 

jU\ 

1-1 \ 

Ul 2 

Uu 

Vl2 

Vll 

\ui 

1-1 

U 2 \ 

u 22 

V21 

v 22 

\u 2 

1-2 

u 22 

U 2 l 

V22 

v 2 \ 

\u 2 

1-2 

0 

0 

0 

0 

Ui-\-Uu~\~U\ 2 

yi+yu+yu 

V° 

0 

0 

0 

U2-\-U 2 i-\-U 2 2 

1 / 2 + 1 / 21 + 1/22 i 


The solution to this problem implies the diagonalization of a 6 x 6 matrix, which is a hard 
problem. Nevertheless, if we decompose U^ and V^ v into irreducible representations 
using Eq. (fT9l) and so the relation X^ u = \rj ,w X a a + X'^ + we would obtain the 
three linear systems 



Uii-\-U\ 2 -\-Ui 
U 2 i + U 22 + U 2 


vn+v 12 +vi V /7° 0 \ 
V 21 +V22+v 2 JyV a a J ’ 


I u n - uu v u ~ v 12 

\u 2 i—U 22 V 2 \ — v 22 


U?\ 

K v y 


Un + Ui2 Vu + V \ 2 
u 2l + U 22 V 21 + v 22 



(11a) 

(lib) 

(11c) 
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which are much easier to solve. The convenience of this procedure is more noticeable 
when we analyze the /3-functions of higher-rank tensor fields. The notation introduced in 
this section will we used in the following ones. 

4 Solutions to the /^-functions of a Yukawa theory 

4.1 Lorentz invariant couplings 

Since the /3-functions of the Lorentz-violating parameters depend on the parameters as¬ 
sociated with the Lorentz invariant theory, we should find the scale dependence of the 
latter first. The /3-functions associated with Lorentz-invariant coefficients do not have 
sign restrictions because their observed values do not vanish at low energy scales. As 
indicated in [30], neglecting 0(c a a , K a a ) contributions, the couplings g and g' are given 
by 


g(P) = gF(p ) 1/2 . g'(p) = g'F(p) 1/2 , 


( 12 ) 


with 


F(p ) = 1 — 5r]g+hip‘ 


(13) 


We will use the definitions r/” 1 = lQn 2 and g± = g 2 ± g' 2 \ the p-dependence of F will 
usually be dropped. In this theory, the values of both g and g' grow as the energy scale 
increases and present the usual Landau pole as F = 0. Fortunately, the Landau pole 


takes place as pl = exp 


|pr , which is a very large momentum scale. For g — g' — |, 

for instance, p^ ~ 10 54 , so the perturbative regime includes a very comprehensive energy 
spectrum. 

Having found the scale dependence of g and g', we can now solve for m and m!. In 


terms of the eigenvalues A± = g 2 + ± y/ (5 g 2 + 3 g' 2 ) 2 + 4 g 2 g' 2 ) we obtain 


87T 2 1 

(d) - ( d ) tp~TK X± — - ( d ) TP 

m± = m± b a + =m±r± 


Hence, in terms of F^> = F + ± F_ 

m - L mF M + m(5g 2 -3^ 2 ) + 8m^ ( _ } 
2 + 4t7 _1 A + — 2g\ 

m > = I_ ™'{hg 2 -Sg? 2 )-mgg' ^ 
2 4r/-U + - 2 g\ 


(14) 

(15a) 

(15b) 


In order to find the momentum dependence of the boson mass and the scalar coupling 
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Sip 2 



Figure 1: Boson mass as a function of p. 

Here pi = JL s- In the solid lines we used 
m 0 

the conditions m = m! = mo, g = 0.2, 
and g' = 0.1. For the dashed lines rh = 
fn! = m 0 , g = 0.1, and g' = 0.25. mo is 
an arbitrary mass scale. 


A/A 



Figure 2: Scalar coupling as a function 
of p. Here An = — • In the solid lines we 
used the conditions g = 0.2 g 0 and g = 
O.lpo- F° r the dashed lines g = 0.1 g 0 and 
g' = 0.25(/o- go is an arbitrary constant. 


A, we need to solve, respectively, the differential equations 

.dp 2 


P- 


dp 

dX 


= 8p 


p g + — 2 (mg + mg 1 ) — ( m 2 + m z )gp 


P~p = V 
dp 


3A 2 + 8A gl - A8g 2 


(16) 

(17) 


where O(c^) and O^K^) have been neglected. Since we have not been able to obtain 
analytical solutions for the last set of equations, some numerical solutions are shown in 
figures Q] and [2j 


4.2 Lorentz-violating parameters 

Now we will solve for the running couplings associated with the Lorentz-violating theory. 
We will independently focus on each one of the four groups presented in table [T] 


4.3 Solution to group 1 

For this family we have the set of differential equations m 

(PkY" = 4 r)g\(K» v + c“ rT - - c^), (18a) 

(&r = 1 sV9 2 +{c^+ - K^+ (c a a +K a a )r)^). (18b) 
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In order to solve previous equations we will use the rank-two irreducible decomposition 
X% = Xg = \rTX a a . ( 19 ) 

Using this decomposition we find that (/3k) a) = (Ye) 1 a = 0; this is expected in the case of 
K /JU because it is symmetric by definition. Writing our linear system as in Eq. (J5J) with 
(; x T )^ u = (K'Y, cf) (and something similar for the traces) we find the solutions 

I& )a a = K^ a a F~ 3/5 , c {d)a a = Y d)a a , (20a) 

K (d)»v = 7/i5 5 c (d)»v = -{d)^ ( 20b ) 

So, in the original basis 

K° a = -\c a a { 1 - F" 3 / 5 ) + \R a a { 1 + 2F- 3 / 5 ), (21a) 

c“« = \c a a (2 + F - 3/5 ) - \K a a { 1 - F~ 3 / 5 ), (21b) 

Kg' = yhf(l - F" 7/15 ) + \KY{1 + 6 F" 7/15 ), (21c) 

cf = (6 + F~ 7/15 ) + yFf (1 - F~ 7/15 ). (21d) 


4.4 Group 2 


Having solved for the coefficients cand K^, we proceed with the vector parameters, 
which are <F, e M , / /J and IF Since all of them are rank-one tensors, they do not need to 
be decomposed into irreducible representations. In this case, however, the /3-functions of 
the linear system do not explicitly depend on cF. Consequently, we must exclude cF from 
the linear system, which is not invertible, and find its scale dependence after determining 
the other three solutions. Thus we must solve the set of differential equations m 

( 6.4 • 

-v 

-2 g 2 gcf 


= V 


-g'gl 

2 g' 2 


-ggl\ ( F J \ 


V 


2 gg' 
2 g 2 


F 


V 7 


e (d)n _ g(d)/x 


(77 
\iP*Yj 

which generate the solutions 

j(d)n _ J{d)ii j{d)fi _ j{d)f! p 

with F±(p) = f _ ( 4=fv ^)/ 10 . 

In the original basis we have, with F^ = F + ± F_ 

JM = + + r /p _ 2F)F(-\ 

r = i-fsr - s'eo + -UL F (-) + g'fp"+g'/") f2F( -) 
’ sl (3 3 1 Vsgl °^’ 2 1 


( 22 ) 




■9 J/ ' F (-) + 9&L_ 


2VQg\ 

g'F) 


76F (+) ), 
(2F ( - } + 76.F (+) ). 


(23) 

(24a) 

(24b) 


276.7 


(24c) 












On the other hand, can be determined by solving [30] 


J i p 


4:(gm+g'm')P l — g^me^+m'f^] 


(25) 


using Eqs. (fT2lh tll5a[l . ( 115bl) . and fl24aD — (124cH . Although a^ can be solved analytically, 
its expression is quite long and so it will not be shown here. An interesting fact though 
is to study its leading contribution at low momentum (p <C 1); it is given by 


a^(p) ~ a M 


+ dF(l-F a ), 


(26) 


where a = ^ ^1 + 2\f§ + (5 g 2 + 3 g' 2 ) 2 + 4g 2 g ,2 j and A is a constant with units of 

mass. The maximum value of a is achieved as g' = 0, thus ct max = ^j(3 + \/6). In spite 
that a > 0 and so increases as p —> 0, we do not really expect to observe large values of 
at low energy scales. For g = 0.25, g' = 0, and p = 10 -100 , for example, F“ max ~ 1.42. 
So = 0(a M , A/ M ) at such low scales. 


4.5 Group 3 

In order to solve for d^, H^ and ( H and U lv are antisymmetric by construction), 

we notice that they obey the set of differential equations m 


(A*r = 
(M7 = 

(Pl)7 = 


\gl( 5d^- d^- V ^d a a )-2g'L^-2gLT 


g 2 +H%'+ 4{gm + g'm')L (f+4 (mg' - m'g)L*^ u 
2(mg 2 _ + 2m! gg')d*jp u + 2(m! g 2 __ — 2mgg')d^ 

W+ 


' 2 7L(f + 2gd*Jp u - 2r/d f 


(27a) 


(27b) 

(27c) 


where we used the Hodge dual representation X*^ v = \e^ val3 X a p. (From now on we will 
use the convention £ 0123 = —£0123 = 1.) After decomposing d 111 ' into its irreducible parts 
we End 


d a a = d a a , d% v = dt! u F~ 2 / 15 . (28a) 

The antisymmetric part of d^ p couples with Eqs. (127bl) and (127c[) and satisfies 

OW = HaK ~ a'L7 - gLT ) • ( 2 ^) 

In spite that Eqs. (127b(127c[) . and (I29a,p only involve antisymmetric representations, 
the Hodge dual components are also present. Such system cannot be solved directly 
because, although the Hodge dual representation of an antisymmetric tensor carries the 
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same information and should have the same scale behavior than the original tensor itself, 
different components are mixed. For example, for n = 0 and v — 1, L^ A = L° A and 
L* A V = L 2 A . This would mix different components of L^, for instance, in Eq. ([29a]). 

The solution to this issue is solved by taking the Hodge dual of Eqs. (I27b| ) . (I27cj) . and 
fl29a,p to obtain a six-variables linear system. In this case we can define the six-component 
vector ( x T Y v = (dn, d* A , H Al H A , L A , where, for example, we have the condition 

(/3d)* A u = ‘^v(g+d* A l/ — g'L* A v + gL A ') from Eq. ([29a]). We thus find the solutions 


d 


H 


A 

(d)nv 


= d^F u Flip) = f 1 / 10 , 


a =H^F 2 , f 2 (p) = f- - , 


9-V2 

20 


l ^ = l ^ f ,, = f 


9+a/2 

20 


Using previous definitions and = F 2 ± F 3 we find 

^(fLT+S’L^), 


JF _ 1 
a A ~ 2 a A 


p(+) I 5F 23 ^ 

^23 “T 


y/41 


HT = H a Fi + ± (2 m-gg'-m'g 2 _)d^-(2m'gg'+m?_)dT 

L 


AF, 


(-) r 


U 4 T g 2 . 


{mg + m'g')L^ A + {mg'- m'g)Lj u 


TfJ-f _ 1 rl 


I?( + ) 5 T23 ) 

,^ 23 vTT 


+ 


2F, 


(-) 


v 7 !! 


(g'dZ-gdT). 


2F, - U 


(+) 

23 


(30a) 

(30b) 

(30c) 


(31a) 

vTT 

(31b) 

(31c) 


4.6 Group 4 

This set of differential equations is the most difficult to handle; it couples the terms 
g XtJ,l/ , J M , and as well as their dual representations. The difficulty mainly arises on 
finding an appropriate irreducible decomposition for g A/ ^; a careful decomposition is per¬ 
formed in appendix [A] An anti-symmetrization in the indices {A,/i} of g Xlw must be 
performed so it can be written according to Eqs. (I67ap — (I67dl) . 

Using the relation = (A,)7 + + \g Xu {0 a ) 2 - and after 

decomposing into irreducible representations we have, in terms of pbV]" = 

m 


(A,) v " = 2 rj {g 2 gT+ gJ’^+lgg'gX - fas'aA" +1(2 g 2 + a' 2 )aA -fa 2 g F 

+ t? AM (,j(a 2 + 6</ 2 )92 + fsiisi + '/d)'‘ + aa'alU 1 ” - laa'a!-** + A'A'' 1 '' 

(32) 

where X a M m = X Al/#1 — X^ vX . Multiplying last equation by eA^o- we obtain 


GOT = v ^ 2 gJ + ^ 


(33) 
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where we used the condition + g% + g% = 0 (see appendix [A] for more details). 

The j3 -function of g%, as defined in Eq. (I63]h can be found by multiplying Eq. (j32]) 
by r)\ v . In order to find (/3 g ) j* we must perform index permutations in Eq. (j32|) before 
contracting the hrst two indices (this must be done because Eq. (j52j) does not provide 
enough information). After doing this we find 


(« = 
(M = 


V+ V 2 )^ + k- 9 i + 6 g'J»+ 6 gg'gj 


(2 g 2 + g' 2 )g f + \glg% - 3 g'J*- 3 gg'g^ 


/ */J, 


(34a) 

(34b) 


The /3-functions of and were already found in [30j. Defining (x T Y = ( 6 , J, gi, g 2 , g\Y 
we obtain, according to Eq. (JHJ) , the matrix 



g\ 2 (mg+m'g 1 ) 

0 

\m'g 2 + 

\mg\ 


0 

8 g\ 

0 

g'g 2 + 

3 gg\ 

V 

0 

-3 g' 

|(2 g 2 +g' 2 ) 

b 2 - 

-3 gg' 


0 

6 g' 

¥- 

i(<7 2 +5 g' 2 ) 

6 gg' 


1° 

2 g 

m 

1 gg' 

2 g 2 


(35) 


While one of the eigenvalues (the one associated with 6 M ) of the previous matrix is always 
Ai = gg 2 , and so is positive, the other four eigenvalues take very complicated forms and 
their signs depend on the values of g and g'\ this behavior is shown in figure|3] For some 
values of such couplings we have a negative eigenvalue and even complex results. The 
complex results associated with A 3 and A 4 in figure [3] are, however, just a mathematical 
artifact and the final results are real in the original basis (this will become more clear 
below); since their real parts are always positive, they do not generate an asymptotically- 
free behavior. 

In spite that A 5 is negative for most values of g and g ', the evolution of the parameters 
described by the matrix in Eq. (135j) also depends from electromagnetic contributions [21]. 
Hence, the inclusion of electromagnetic effects may change such tendency. Additionally, 
the magnitude of A 5 does not take considerably large values for typical values of /3 (this 
parameter is defined in figure [3]). In figure [3] the smallest value of A 5 is A 5 ~ — 8.8 x 1CT 3 
and takes place at /3 = 3. Using this value and g = 0.25, we see that F~I A5 I ~ 0.96 for 
p = 10~ 100 , which is very close to 1. Therefore, A 5 induces a negligible deviation at low 
energy scales. 

If we neglect electromagnetic effects, there are still some conditions that could make all 
eigenvalues in Eq. 1)35]) positive for arbitrary values of g and g'. The solution is assuming 
that the values of some coefficients at the renormalization scale vanish, so there is no 
LV at any energy scale coming from such fields. Nevertheless, such conditions are very 
unlikely and restrictive and demand fine-tuning. 
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Figure 3: Values of the eigenvalues A 2 , A 3 , A 4 and A 5 as a function of g and g'. Here g = g 0 
and g' = /3g 0 - The vertical dashed lines divide a set of intervals; the horizontal dashed 
line is the value 0. The first interval is 0 < (3 < and all eigenvalues except A 5 are 
positive, like in the third interval given by 1 < /3 < 1.038 and the fifth, where j3 > 1.434. 
In the second interval all eigenvalues are positive and < /3 < 1. In the fourth interval, 
1.038 < /3 < 1.433, A 3 and A 4 are complex. The upper small graph zooms the interval 
where A 5 > 0 , the lower small graph zooms the imaginary parts of A 3 and A 4 . 


Using Eq. (|35li into Eq. fl32l) and the /3-function of g Xui1 (last one can be found by 
permuting indices) we have that the other two irreducible representations satisfy 


(A,)U = b[(2 9 2 +s' 2 ) 9l A r- 9 ! 9 A r - 9s'j!- i *+ 2 99 ' 9 /y-' 

(ft) 2 A r =| 9 (2 9 2 + 9 ' 2 ) 9 A r- 9 2 9 U - 9 s'j!' a *+ 2 99 ' 9 / + a ‘‘ 1 "' 


(36a) 

(36b) 


Last system is, however, very difficult to solve. For convenience we will change to a differ¬ 
ent basis by defining g± ll/ = g x ^ v ± g\+ v ■ In terms of ( x T ) Xfll/ = (g^, g+ J ' u , g < ^ x ^ l/ 

and the matrix 

(3g 2 0 0 0 \ 

0 3g 2 0 0 

0 —3gg' g 2 + 2g' 2 


M = 


2g 


3 gg’ 0 


-.9.9 

2 


99 

g 2 + 2 9 


1 2 


(37) 




we find the eigenvalues A !)2 = 2gg 2 and A 3i4 = g 2 + 2g'~ =|= igg'. Although we have two 
complex eigenvalues, the evolution of the given parameters is real when we go back to the 


12 





























original basis (this also applies for the parameters described in figure [3]). In the diagonal 
basis we find 


9+ 


_ -(d )\fiu T-, 

— 9- *i, 

_ -(d)Xfiu rp 

- 9 + h 2 , 


(d)[*X/j]v _ (d)[*\/i]u jp* 

9- ~ 9- r l ) 

9+ ~9+ *2 > 


Fi(p) = F~~ 92/5 ~ 9 + , 

F 2 (p) = j p-(9 2 + 2 ff' 2 )/( 15 4) e xpfi&^ 


(38a) 

(38b) 

(38c) 

(38d) 


In terms of F 2 R = Re [F 2 ] and F 2 = Im [F 2 ] we have 


A/xi/ _ 1 -\fiu jp , 

9\— - 2 9- + 2^+ 


1 vA/XZ' TpR 


1 Tpl I 

~o9+ b 2 + 


3_ 2 _, 2 _M, + 6 j 2 


8(jP 2 + 2g 2 g' 2 


F 1 

r 2 


0 _2_/ 2 _Aui/ /->- 2 —/_l*A/i]y 

3 g g g- - 6 gigg'gL 


Xiw A av 

g 2 + =9 1- 


8(^)2+ 2^/2 
■ S-^Fi. 




(39a) 

(39b) 


At first glance, one might be tempted to think that and g^+ u have an oscillatory 
behavior because of the definition of F 2 in Eq. (138d|) . Nevertheless, such oscillatory 
behavior does not take place. From Eq. (fl3|) we can easily see that | In F| does not take 
very large values unless p is extremely small or large. For instance, for g — g' = 0.25, 
lnF(10 -50 ) ~ 0.65 and lnF(10 50 ) ~ —2.42. On the other hand, if g' = / 3g , ~^r is 

maximum at (3 = 1, hence ^||t < A. Therefore, we can easily conclude that |^|| T lnF| < 
| for typical values of g, g', and p. 


5 Application to Lorentz-violating QED 

The method that we have used to decomposing into irreducible representations can also 
be applied to Lorentz-violating QED using the results from [3T] . 

In Lorentz invariant QED the evolution of the electron charge and its mass are given, 
respectively, by mm 

e 2 = e 2 Q~\ m = rhQ 9/ \ with (40) 

Q(p) = 1 — (e 2 /67r 2 ) lnp. (41) 

From pH] we can see that the parameters describing Lorentz-violating QED can also be 
decomposed into four groups, which are 

1. Set {g^c^, {k F ) x ^ a }. 
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Figure 4: g X ^ u as a function of the scale. 
Here the black lines represent the condi¬ 
tions 2g = g' = 0.2g o and the gray lines 
the conditions g = 2g' = OAg^. In the 
continuous lines g X ^ u = 2g X ^ u = 2 H 0 ; in 
the dashed lines g X,w = g^ v = A 0 . 


Figure 5: g^+ u as a function of the scale. 
We used the same initial conditions and 
labeling as in figure |4l 


2. Set {a M , e M , / M }. 

3. Set {dJ LV ,W lv }. 

4. Set {b x g llu : g Xflu ,(kAF) x V^ u }- 


5.1 Group 1 

In Lorentz-violating QED we also have the parameter ( kp)\^ a associated with the elec¬ 
tromagnetic held. Since its contribution to the Lorentz-violating QED lagrange density 
is A C = -\{k F )x^v ^F x > 1 F ua , it has the symmetry properties [3T] 


{J^f)X fiva ( kp')vcrX/i {kp)^,\ua {kp'j Xficru i 

[kp) XfM/a "F (kp)xuafi 4“ (^If) Xcrfiv 0 , 

(k F ),r = o. 

From (Mj we have the /3-functions 


(A 


kp 


\X[iva _ 


67 r 2 
*2 


(k F ) x ^ a - r} Xu c^ ] + g [ ^cg a] - \c a a g [Xv g^ ] 


(pcT=h 2c r + W a c a a - c k F r v 


(42a) 

(42b) 

(42c) 


(43a) 

(43b) 


where we used the definition X^ Xu Y^ = x Xu Y^ a —X^ U Y X(T . Notice that previous relations 
satisfy the symmetry conditions stated in Eqs. (I42al) — (I42cl) : these symmetry conditions 
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also imply from Eq. (143bj) that c^ u = c^ u 


From now on we will use the same definitions 
as in Eq. (1691) and deal with ( kp) x ^ ua as explained in appendix [B] Notice that when we 
decompose {kp) Xtlua according to Eq. (1731) we have permutations in the indices {A, q, u, a} 
just like in Eq. (143 al) . Hence, each permutation of indices satisfies the same equations 
and so a decomposition of the /3-function into irreducible representations can be written 


as 


(#,) 




= h 


VW.+MOf = A 


(k F ) x r° + V" 


«Mr 


2 cn 


loiv 1 


2cr+5<r<=“ - (ms 


(44a) 

(44b) 


Taking the double trace of Eq. (I44a|) we see that c a a = 0. This is not inconsistent, 
however, because this parameter has no physical relevance because it can be absorbed by 
a redefinition in the fermion held |38j . Multiplying Eq. (I44a,p by r)\ v we have the linear 
system 


0 8k F )%\_£_( 1 -2 

mr) \-l 2 J{(%' 

The solutions are thus 

(MS = (MS - 1 ((MS - 244 (i - Q ~ 3 ). 
4’ = cr+K(Ms-2cS’)(i-eN' 

Previous results imply that 

IMP" = (MP'"Q _1 • 



(45) 


(46a) 

(46b) 


(47a) 


Notice that the solutions found in Eqs. 046al - l4Ta]) are the same than those found in [31] 
once we write back each tensor held in terms of their irreducible representations. 


5.2 Group 2 

As stated, vector parameters do not need to be decomposed into irreducible represen¬ 
tations, so the decoupling is direct. From UD we can easily see that e tJ ‘(p) = and 
/ M (p) = f 11 . We use such results to integrate a^; notice that this system cannot be solved 
through a matrix diagonalization because (/3 e ) M = (/3/)^ = 0 (the matrix would have two 
rows of zeroes). We confirmed the result [34) 

a" = a^-me^l - Q 9/4 ) . (48) 
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5.3 Group 3 

In order to find the solutions for this group we first decouple d liU , whose (3- function does 
not depend on H^ EH. The solutions are 

d a a = d a a Q~\ d^ = d% v Q-\ d% = d%. (49) 

For H ,,v = H 1 ^ we use the solution d*/ v = d*^ v to integrate H'Y . Like in the previous 
group, we have to follow this method instead of the matrix-diagonalization-procedure 
because {fid) 1 a — 0; the integration gives 

Hfifi = H^Q~ 3/a + rhd*/ u Q 9/4 {l - Q~ 3 ). (50) 

While the solution for d^ is the same, the solution for Hfifi differs from the one found in 
[51] : the reason is probably a typographical or algebraic mistake. 


5.4 Group 4 

We will first focus on the solutions for the coefficients g^ [31]. By permuting indices and 
decomposing into irreducible representations we find the /3-functions for g x ^ v and g Xufl , 
which are, respectively 


wT +!(&)£“• + yy/y; = 
-u>,jT+= 


e 2 
87 r 2 

e 2 

8tt 2 


f 91 -" -1 dr++ W 
| 92 A r+ la'r+ +/y M 9 f 


(51) 

(52) 


where we used again the condition g^ + g% + g% =0 (remember X A ^ ,i = X Xu ^ — X^ vX ). 
If we multiply by £\ lwa in either equation we see that ( fif )* M = 0. Thus g X fi w = g X fi w or 
g*Jfi = g*fi. Multiplying by rj\ u in both equations we find the linear system 


(fi 9 Ti\ 

(« J 




(53) 


leading to the solutions 

= 9i‘(5<2- 8 '' 8 - 4Q- 3 ' 2 ) + 2g»(Q-W - Q ~ 3/2 ), 

= 9?(5Q- 3/2 - 4Q-W)+ 10jf(Q- 322 - Q- s/s ). 

The scale dependence for the other rank-one coefficients is given by 


(k AF r = {k AF YQ-\ 

b' 1 = b»- mg*f( 1 - Q 9/A )-\{k AF y{ 1 - Q- 1 ). 


(54) 
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(55a) 

(55b) 







lOW/go 



Figure 6: A 6 as a function of g , g', and e. Here we have g = g 0 , g' = (3go, and e = ag 0 . 
For each one of the six plots we have a\ = 0.25, a 2 = 0.6, cr 3 = 1, cr 4 = 2.5, a 5 = 6, and 
ctq = 20. It can easily be seen that when electromagnetic effects are included in Eq. ([35]) . 
the negative eigenvalue tends to become more positive. This effect is larger for larger a. 


Additionally, using previous results we obtain the linear system 


(ftd-T = jLf 4 -A f AT 
our; stt 2 y i 


(56) 


whose solutions are 


gX,u = g^ 2Q - 9/4 _ Q - 3/2) _ 2 g*r{Q~ 9/4 - Q~ 3/2 ) , ( 57 a) 

g x 2 r = g x r{2Q- 3/2 - g- 9 / 4 ) - g^(Q- 3/2 - g- 9/4 ) • (57b) 


5.4.1 Inclusion of Yukawa terms in g x ^ 

Let us remember that Eq. (155[) described the set of differential equations for the parame¬ 
ters J^,g^,g 2 , and g*/ in the pure Yukawa regime. When we also include electromag¬ 
netic contributions we will find a linear system of the form ((3 x )m = El= i M mn x%, with 
[x T Y = (k A F, b, J, gi, g 2 , g\Y and 



(V 

0 

0 

0 

0 

0 

\ 


6e 2 

g 2 + 2(mg+m'g') 

0 

\m'g 2 + 

| mg+ — 6 me 

2 

V 

0 

0 

00 

+ 10 

0 

g'gl 

3.9.9+ 


0 

0 

-3 g' 

\{2g 2 +g' 2 ) — e 2 

1 9 2 - - 2e 2 

-3 gg' 



0 

0 

6 g' 

\g 2 _ + 10e 2 

l(g 2 +5g /2 )+8e' 

6 gg' 



(° 

0 

2 g 

1 99' 

1 99' 

2 g 2 

/ 


Last matrix has 6 different eigenvalues. While two or them can be read automatically from 
last matrix and are Ai = ^e 2 and A 2 = gg\, the remaining four take very complicated 
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forms. Nonetheless, three of them, name them A 3 ,A 4 and A 5 are always positive for 
any combination of g, g\ and e. Xq is not always positive for any combination of such 
parameters, which can be seen in figure [ 6 ] 

6 Summary and Conclusions 

We have used a method based on decomposing a tensor into its irreducible representations 
to solve the /3-functions of two theories with Lorentz and CPT violations. Although we 
mainly focused on a Yukawa theory because the scale dependence of its coefficients had 
not been solved, we also applied this method to an electromagnetic theory for comparison 
reasons. Indeed, up to some possible algebraic or typographical mistakes we confirmed 
the results found in [3Tj . Nonetheless, our method is more convenient because it reduces 
the complexity of the problem. 

As stated, those operators whose scale dependence is described by negative /3-functions 
could lead to an inconsistent theory if their scale dependence is strong enough, thus 
predicting a large imprint of Lorentz violation at low energies, which, at best of our 
knowledge, has not been observed. Fortunately, in their respective diagonal bases, most 
coefficients are described by positive /3-functions. In this way we could envision, within 
the scope of these models, Lorentz symmetry as an emergent low-energy property. 

Some mass operators have negative /3-functions; however, their scale dependence is 
weak enough and so their imprint of LV at low energies is negligible. This fact also took 
place in some of the irreducible representations of g XfJ,u . Since the solutions obtained for 
the analyzed Yukawa theory do not predict large imprints of Lorentz violation at low 
enough energies, this model is consistent with observations, at least, at one-loop order. 
Actually, the values of the couplings and masses of the Lorentz invariant theory are not 
constrained by this condition. 

This method could also be applied to renormalize other Lorentz-violating theories with 
similar characteristics or any other theory that has similar properties. Studying the scale 
behavior of a Lorentz-violating QCD theory could be of particular interest because, in 
the Lorentz invariant case, it is asymptotically free. 


Appendices 

A Irreducible decomposition of three-rank tensors 

The permutations of three indices form the group S 3 . The irreducible representations of 
three-rank tensors can be found using the following set of Young diagrams 


18 


s 


A 


1 ~ 


V 


A 

A 

V 


1 + 


A 

f* 

V 



A 

V 




Last set of diagrams generates the four irreducible representations 


X x ^ v = ~(^X ( ' X ^ >1 ' + X x ^ tl + X u ^ x ^ 

(59a) 

X^ v — I 

(59b) 

X x £ v = i(X ( - A/ ^ 1 ' — X v ^ x ^ 

(59c) 

X Xfll/ = h(^x' [X ^ u + x^ A ) 

(59d) 


with X^ u = X Xtll/ + X ^ Xv , = X x ^ v — X mAi/ , X a ^ = X Xufl + X fluX , and X x ^ = 

X Xufl — X^ uX . However, last set only generates 4 linearly-independent representations. In 
order to complete six independent combinations, we will exchange the indices /j, — > v in 
the tensors generated by last two diagrams. So we define 


-\rXf-LV _ -yXlSfl -yXfllS _ -yXufl 

A 2+ — A 1+ , A 2 _ — A x _ 

This set of six tensors can be inverted to find the relations 


/ X x > w \ 


f 3 3 4 0 4 o\ 


(X^\ 

X Xv v 


3 -3 0 4 0 4 


\rX^LU 

X^ vX 

1 

3 3 -4 -4 0 4 


J^X [iv 

X^ Xv 

X vX v 

^X^J 

“ 3 

3 -3 4 0 -4 -4 

3 3 0 4 -4 -4 

\3 —3 —4 —4 4 O) 


\rXfliy 

a 2+ 

-^X^LV 

Ui""/ 


We can verify that X* Xfll/ = X*^ u as well as the following relations 

t /" flXlS ~^rX[AlS ^rXl/fl ^riyXfl \^X f^lS 

A 1 + = A 1+ , A 1+ — A 1+ = A 2 _|_ , 

- (x A r+ w 2 a d, 

^rUflX ^rXufl ^rfllsX ^rX^LU 

A^_ Ai- ? Ai_ A^_ A 2 _ , 

H’Xv X^Xl 1 ^j^Xiiv j^-Xiii'^ 


■\rVfiX _ \rfllSX _ 

^1+ — ^1+ _ 


(60) 


(61) 


(62a) 

(62b) 


Now we want to isolate the “trace-like” components, which have different transformation 
rules. Let us define 


XI = 1° 


x° 2 = X Q " a 


XI = X c 


(63) 
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In order to isolate these contributions we must make the six irreducible tensors previously 
defined traceless in all of their components. So we make 

X*r = - X r } Xl \X x + X 2 +X 3 y - X V x ^ (X 1 +X 2 + X 3 y, (64a) 

X A f = X x /\ (64b) 

x x y = x^r+yyix, - x^y-y^yx, - x 3 y, (64 C ) 

x x y = x x y v +yy\x x - x 2 y-y v x ^\x 1 - x 2 y, (64d) 

where the subindex “t” means that the tensors are traceless and, again, X Xl X = X x ”+ 
and X x y = X x y. The traceless representations of X x ^ u (and its permutations) can be 
found using Eq. (T6TT) . 

From now on we will drop the traceless notation “t” and assume that the irreducible 
tensors in Eqs. (I64a,p — (I64d|) are traceless. 

For our particular situation of interest, we will focus on the properties of X XflL ' when 
it is antisymmetric under the exchange A <<->■ /i (this is the symmetry condition obeyed by 
g x fiu ag i n t roc i ucec i i n Eq. fj2a|l ). This means that we can write X XfJ,u as X XfJ,u = 

This condition reduces the number of linearly-independent irreducible representations 
(neglecting the “traces”) from six to three, which now are 

X [Xtt]u = | ( X [Xfl]u + X^ x + X [uX] y , (65a) 

x \Mv = 1 - l -X [tlX]u + \x u[x ^ + \y [v] {X 1 -X 3 ) M ), (65b) 

X [x Y ]v = \ (X [Xfl]u - \x^ x + ix u[flX] + ( x ± - X 2 Y) ■ (65c) 


Notice that V\„X [X ^ = 0, r jXu X^ = yX 2 -X^, and ■q /u/ X^ v = |(X,-X 2 ) A . We will 
now use the definitions X-^ = |(X“ — X 3 ), X^ = |(A"^ — X"), and xY = |(X 3 — X%). 
By construction we have the property 

x [ y + x [ y 1 + x [ y = o. (66) 


In order to simplify the notation we will assume from now on that X Xfll/ and its permu¬ 
tations are antisymmetric under the exchange A -H- /i and so X^ 1 ' —>■ X Xflu . Relabeling 
xf 1 x ? we find 


X x ^ = X H*V + 4 x a ^ + LjjWxg, 

x Xuti = -x x y v + \x x y v - l -rj x[v] x>\ 

X^ x = -X x r - ix A r + \x x y u - \r) X[v] X>Y 

j^/iuX _ _ i'll _ _ -^X[iv y^vX[i _ _ -^ufiX 


(67a) 

(67b) 

(67c) 

(67d) 


where X" satisfy Eq. (J66]). 
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B Irreducible decomposition of four-rank tensors 

A four-rank tensor X Xp,va must also be decomposed into irreducible representations under 
the group S 4 ; its associated Young diagrams are 


2 

P 

V 

o 


2 a 


2 

P 

V 

o 


2b 


la 

2 

P 

V 

o 



2 

V 

P 

o 


3 a 


2 


v 


a 


P 


lb 


2 

p 

<7 

V 



3b 


2 


P 


o 


V 


lc 



2 

V 

a 


P 


3 c 

A 

2 

o 


2 

P 


P 

V 


V 



o 


and permutations of the indices {A,/i, v, a} to complete 24 different representations. In 
terms of the set of Young diagrams we can write 


X xpya = x 1 X XtxVfT +x 2 X x r a + x 3 X x r+ x 4 X Xpva + x 3 X xpua + x 6 X Xpuo 

+ X7 xir + Xs x x r + X9 x x r + 


( 68 ) 


where Xi are rational numbers. Fortunately, the /3-functions of Lorentz violating QED 
that deal with four-rank tensors only include the combinations (/ cf ) Xp.ua and (/ cf )\ ap a - 
Hence, we do not need to End all irreducible representations, we just need to focus on the 
trace and double trace representations of (kp)\^ V(T . Let us define 


afU 0 _ ^pa f) — l£P a 0 


•spafi _ vapf) 

A 4 — A p 


x 7 =xr 


p p' 


xf = x ap / 

A« = X pp ' 


pp 


\ra/3 _ v&Pp 

A 6 — A pi 

Xg = X PP ', . 

P P 


(69) 
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In terms of these definitions we can write the traceless and double traceless tensor x x>J ' l ' cr 
as 


fiver _ ^X^iua 


96 ' 

— n Xv 
— n X(T 

96 ' 


- 31Xi+X( +7X 2 —X 2 -X :i + 7X‘-X A +7Xi + 7X 5 -X 7 -3X 6 -3X 6 ' 

1 M* 7 

7X 1 -X 1 r -31X 2 +X 2 T +7X 3 -X 3 T -3X4-3X 4 T + 7X 5 -X 5 t -X 6 + 7X 6 t 

- Ad + 7Xf+ 7X 2 - X 2 T - 31X 3 + Xl + 7 X 4 - Xl - 3X 5 - 3X 5 T - X 6 + 7X\ 


) 1 U 


96 

1 

96 

1 

96 


V' 


p,is 


7Xi - X[ + 7X 2 - X 2 J - 3 X 3 - 3 X 3 + 7 X 4 - Xi - 3 1 X 5 +X 5 r +7X 6 - X, 


Act 


n 


v 


1±<j 


- Ad + 7Xf — 3X 2 — 3X 2 + 7Ad - X[ - 31X 4 +Xj + 7X 5 - X 5 T + 7X 6 - X 7 

- 3X x - 3Xl - X 2 + 7Xl - X 3 + 7Xl +7Ad - Xl + 7 X 5 - Xl - 31X 6 +X 6 7 


1 Xu 


Xfi 


—ri x ^ri ua 
144 ' ' 

—v Xa v^ 

144 ' ' 


17X7 - 7Xs - 7Xq 


+ i^ A V" 


- 7Ad + 17X 8 - 7Xg 


— 7 X 7 - 7X 8 + 17X 9 


(70) 


where (X T ) ap = X^ Q . 

We will now study the properties of four-rank tensors obeying Eqs. fl42a)) — (j42c[i . 
When we multiply by r] X r ] IJL(J , and Whiter respectively in Eq. f!42bl) we obtain the 
conditions 


(Xi)f + + (Ad= 0 , (71a) 

(X 4 )^ + (Ad)^ + (Ad) 5 ^ = 0 , (71b) 

X 7 + X 8 + X 9 = 0 . (71c) 

Eq. (142cl) implies X 8 = 0, so X 7 = —X 9 . However, X 7 = X Xflua rjx^rj^ = 0 because of the 
antisymmetry properties in Eq. (142al) . so X 9 also vanishes. This implies that X x ^ a is 
double traceless in all components. 

Using the symmetry properties stated in Eq. (I42a,p we find that X x ^ va can be simplified 
to 


Y A„,„ = x ^ u , + 1 - V'(xX + ijf’py- 


(72) 


where X a = X 2 - X 3 + X 4 - X 5 . 

Notice that the symmetry conditions obeyed by X x ^ va also imply that (X 2 )^ cr = 
(x 4 )f = -(X 3 )f = —(X 5 )g a = \(X a Ys 1 so we finally have 

X x ^ = xY ua + |h Al/ (A 2 )f - |^(X 2 )^ - (X 2 )f + f iT(X 2 )^ , (73) 


which is enough to decompose (kp)x^cr- 
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